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We present a second-order Green's-function theory of the one- and two-dimensional S = 1/2 fer- 
romagnet in a magnetic field based on a decoupling of three-spin operator products, where vertex 
parameters are introduced and determined by exact relations. The transverse and longitudinal spin 
correlation functions and thermodynamic properties (magnetization, isothermal magnetic suscepti- 
bility, specific heat) are calculated self-consistently at arbitrary temperatures and fields. In addition, 
exact diagonalizations on finite lattices and, in the one-dimensional case, exact calculations by the 
Bethe-ansatz method for the quantum transfer matrix are performed. A good agreement of the 
Green's-function theory with the exact data, with recent quantum Monte Garlo results, and with 
the spin polarization oi a u — 1 quantum Hall ferromagnet is obtained. The field dependences of 
the position and height of the maximum in the temperature dependence of the susceptibility are 
found to fit well to power laws, which are critically analyzed in relation to the recently discussed 
behavior in Landau's theory. As revealed by the spin correlation functions and the specific heat at 
low fields, our theory provides an improved description of magnetic short-range order as compared 
with the random phase approximation. In one dimension and at very low fields, two maxima in 
the temperature dependence of the specific heat are found. The Bethe-ansatz data for the field 
dependences of the position and height of the low-temperature maximum are described by power 
laws. At higher fields in one and two dimensions, the temperature of the specific heat maximum 
linearly increases with the field. 



> 

(N 

O 
en 

o 

o 

-(— • 



o 
o 



X 



I. INTRODUCTION 

In the theory of low-dimensional magnetism the es- 
sential role of quantum and thermal fluctuations, es- 
pecially in the description of magnetic short-range or- 
der (SRO) at arbitrary temperatures, is of basic inter- 
est. Whereas for Heisenberg antiferromagnets the in- 
terplay of low dimensionality and quantum fluctuations 
is important already at T = 0, in ferromagnets quan- 
tum fluctuations occur at non-zero temperatures only. 
The study of low-dimensional quantum ferromagnets in 
a magnetic field was motivated by the progress in the syn- 
thesis of new materials, such as the v — 1 quantum Hall 
ferromagnets ji which may be described by an effective 
two-dimensional (2D) S — 1/2 Heisenberg model)2ii4 the 
quasi-2D ferromagnetic insulators A2CUF4 (A=K, Cs)^^ 
La2BaCu05,''^ and Rb2CrCl4fi the quasi-lD organic fer- 
romagnet p-NPNN (Ci3Hi6N304)^ and the quasi-lD 
copper salt TMCuC[(CH3)4NCuCl3].^ 

The 2D 5 = 1/2 ferromagnet in a field was investigated 
by Green's-function decouplings of first order, -'^^ i.e., by 
the random phase approximation (RPA)-'^^ and the Callen 
decouplingjiS by Schwinger boson theories^Si^ and by 
quantum Monte Carlo (QMC) simulationsi^ Thereby, 
the magnetization^iSiiiiS and the spin lattice relaxation 
rate^'* were calculated. The ID ferromagnet was studied 
by the Bethe-ansatz method, where some exact data for 



the zero-field magnetic susceptibility and specific healii^ 
as well as for the magnetization and correlation length 
were given. '^^ Recently, in the ID model a power law for 
the shift of the temperature of the susceptibility maxi- 
mum with the field was reported and argued from Lan- 
dau's theory to appear in the 2D model, too.^^ Therefore, 
a detailed analysis of the thermodynamic quantities of 
the ID and 2D Heisenberg ferromagnets as functions of 
temperature and field is of interest, also for comparison 
with experiments. 

We consider the 5 = 1/2 Heisenberg model 



H 






(1) 



[{i,j) denote nearest-neighbor (NN) sites; throughout we 
set J = 1] along a chain and on a square lattice. To pro- 
vide an improved description of SRO and of the thermo- 
dynamics (magnetization, magnetic susceptibility, spe- 
cific heat) as compared with the standard approachesjiS 
we go one step beyond the first-order Green's-function 
decouplings. To this end, we adapt the Green's-function 
projection method dealing with se cond t ime derivatives 
of spin operators outlined in Refs. Il6ll7l Furthermore, 
we perform exact finite-lattice diagonalizations (ED) on 
an iV = 16 chain and an iV = 4 x 4 square lattice using 
periodic boundary conditions. 



The exact Bethe-ansatz results for the ID case are ob- 
tained from an eigenvalue analysis of the quantum trans- 
fer matrix of the Heisenberg chain, a concept that is also 
the basis of the work reported in Ref. [l^- Here, unlike 
the treatment in Ref. il4l we perform this analysis by 
solving a certain set of nonlinear integral equations to be 
found for instance in Ref. [3 These integral equations 
were analysed in the literature extensively for the anti- 
ferromagnetic Heisenberg chain. The ferromagnetic case 
satisfies the same set of equations with just a sign change 
in the temperature dependent term. Despite this rather 
minor change in the analytical formulation the numeri- 
cal treatment of these equations is rather different from 
the antiferromagnetic case. The iterative treatment is 
plagued by slow convergence, in particular for low fields 
and low temperatures. A numerically much better con- 
ditioned formulation can be derived by combining the 
methods of Refs. 18 and [I^ Details of these calcula- 
tions will be given elsewhere. Our results are in perfect 
agreement with those of Refs. [ij and[lj if available. 



II. GREEN'S-FUNCTION THEORY 

To calculate the transverse and longitudinal spin corre- 
lation functions we determine the two-time retarded com- 
mutator Green's functions Gq^(w) = {{Sq;S^q))uj {i^f^ = 
H — , zz) by the projection method, where we neglect the 
self-energy»iSJ^ Taking into account the breaking of spin- 
rotational symmetry by the magnetic field we choose, as 
for the XXZ model, ^^ the two-operator basis (S'+,z5+) 

and {Sg,iSq). To approximate the time evolution of the 
spin operators —S^ and —S^, we take the site represen- 
tation and decouple the products of three spin operators 
in —S^ and —S^ along NN sequences {i,j, I) introducing 
vertex parameters a"'^ in the spirit of the scheme pro- 
posed in Refs. [13 andl20l 



S+S+S[- = a^ 



-{S+S^)S++a+-{S+S^)S 



Dj Sj 6; 



'{Sj Si 



>Sf 



(2) 



(3) 



Here, following the investigation of the ferromagnet at 
h — 0^ the dependence on the relative site positions 
of the vertex parameters (cf. Ref. Hal is neglected. We 
obtain 



s: 



[K 



h^]S+ + 2hiS^ 



(4) 



^1 = «)'5^ 



with 



«)^ = -(1 - 7,){A- + 2za-C+-(l - 7q)} , 



2 

(8) 



A- = 1 + 2a-{(z - 2)C+- + C+- -{z + 1)C+-} , (9) 

where Cnm = 2^nm + ^mm ^mn = ^R = i^O ^W ^ 

R = nex + mey, 7^ = f X]i^=i '^^sg^, and z is the co- 
ordination number. This approximation, Eqs. (01) and 
(O , is equivalent to the equation of motion decoupling in 
second order. ^'^ Finally, we obtain 



G+-H = 



A, 



qi 



1=1,2 



ji,2 = h±iu+-, (10) 



■^qi 



IvF„ 

G-H- ^ 



Cj2 - (t^2^)2 



(11) 



with 



^qi.2 = {S^) ± TrT^{M+- - 2h{S^)) . (12) 



2cj, 



The first spectral moments A/+ = {[iS'^ ,S_^) and 
^^q^ — ([*'S'qi S^^) are given by the exact expressions 

M+- = 2zGio(l - iq) + 2h{S') , (13) 



M|^ 



zCt^-{l-lq) 



(14) 



The spin correlators are calculated as C'^ = 
^T^qC^^e"!^ with C^" = (S'^51q). By Eqs. ^ and 
(inj we get 



i=l,2 



f^zz q 

'^ ~ 2w- 



Cl^^Cl^+Dl\ (15) 



[l + 2n«)] 



(16) 



where n{uj) — {e^^^ — 1)^^- As shown in Ref. |22, for 
the complete determination of correlation functions cal- 
culated from commutator Green's functions one has to 
take into account an additional term, if the correspond- 
ing anticommutator Green's function has a pole at a; = 0. 
Here, we haveSi 



^r-ii-|G^^^^^o 



2 



w 



(17) 



The equation of motion for the anticommutator Green's 
(5) function G[^^"{iv) yields Eq. ^ with Af^^ replaced by 



M. 



(+)^ 



2u;G^^ where M, 



A+)^ 



([i5^,5iq]+). Bythe 
spectral theorem for G^^ it can be easily verified that 



i^q ) =;j(l-7q){A+ +2za+ Gio(l-7q)}, (6) a4+)^^ = 0. Thus, Eq. IHIl with tj^^o == yields 



A+- = 1 -H 2a+-{(z - 2)Gii + G20 ~ {z + l)Gio} ; (7) 



^q^ — C'q^<5q,0 — / , C'fl ^q,0 



(18) 
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From Eqs. ((TB|l and ifTH)) we have C, 
relation 



zz 
q=0 



1 9(5^) _ 1/1^ 2 



0. By the 



(19) 



foUowing from the first and second derivatives of the par- 
tition function with respect to /i, in the thermodynamic 
hmit wc finaUy obtain 






iqR 



(5^ 



(20) 



A = hm -—-T, by Eq. (El) with cothtj„/2r = 2T/w„, 
A is calculated from 



iV ^ 1 - 7o 



iV ^ 1 - 7q + A 4T 



(25) 



The zero-field susceptibility is given by Xsi^ ~ *-*) 
{zX)-\ 



III. RESULTS AND DISCUSSION 



Note that the transverse correlator has no additional 
term, i.e. D~^ = 0, because of Wq^o;i,2 = hy^O. 

The magnetization per site, m = —2^b{S^), is calcu- 
lated as 



(S-^ 



Cn 



(21) 



From (S^) the isothermal magnetic susceptibility x = 
4:fi%Xs with Xs = d{S^)/dh may be derived. 

To complete our scheme, the two vertex parameters 
a^^{T, h) have to be determined. To this end, we use the 
sum rule Cg^ — j and the exact representation of the in- 
ternal energy per site u = {H)/N in terms of Gg ~(w),^^ 
i.e. 



z h 1 -s—^ 1'°° cLu 



(22) 



27r 



where e„ 



+ w)ImG'+"(w)n(w) 



Thus, we have 



1 ^ 1(1 - 7q) + h. 
a closed system of equations for nine quantities 
((5^),Cfo^Cl7,C2^o'',a^^a+-) to be determined self- 
consistently as functions of temperature and field. As 
may be easily seen, at T = the exact results (S"^) — |, 
\, C^ = are reproduced. 

In the case h = Q the spin-rotational symmetry, im- 
plying {S^) = and C^~ = 2C^^, is preserved by our 
scheme with a"' = a^^ = a, and the theory reduces to 
that ofRefs.ll^andlli 

It is of interest to compare our results with the RPAii 
which employs the decoupling iS^ — tOqS^ and yields 



C^^ 



Gj-(. 



UJ 



2{S' 



o^^z{S'){l-j,) + h, (23) 



(5^ 



N 



E^oth 



2T 



(24) 



zz 
10 



Note that the longitudinal correlation functions cannot 
be obtained by such a simple decoupling, except for C 
which may be calculated in RPA by Eqs. l|^ and (|^ 

At /i = 0, in Ref. ^^ the RPA was extended to the 
disordered phase, i.e. to T > for ID und 2D ferromag- 
nets (Mermin- Wagner theorem). Introducing the ratio 



A. Magnetization 

Considering the magnetization of the chain, in Fig.^ 
the analytical and ED results as well as our Bethe-ansatz 
solution are plotted and compared with the RPA results. 
Let us emphasize the excellent agreement of our the- 
ory for the chain with the ED and Bethe-ansatz data 
over the whole temperature and field regions. For the 
ID ferromagnet the RPA turns out to be a remarkably 
good approximation for (5^). In the inset the magne- 
tization at low fields is depicted, since the low-field be- 
havior of the specific heat turned out to be of particular 
interest (see below). Note that the experimental accessi- 
bility to the magnetic field strengths B corresponding 
to a given h value may be checked from the relation 
h = 0.116B[T]/J[meV]. Considering, e.g., the quasi- 
ID ferromagnet TMCuC with J = 2.6meV,^ the value 
h — 0.05 corresponds to the magnetic field B ~ IT. In 
Fig. ^ our result for the 2D ferromagnet, together with 
the QMC data for a 32 x 32 system, are shown. Compar- 
ing the ED with the QMC results, the finite-size effects 
are seen to be largest for low fields and at intermediate 
temperatures (cf. Ref. 0); for large fields {h > 0.4) they 
become small. Furthermore, as can be seen in Fig.^ the 
finite-size effects in the 2D model are more pronounced 
than in the ID model. This may be due to the smaller 
linear extension of the 2D system as compared with the 
chain of an equal number of spins. In two dimensions, at 
low fields {h < 0.3) the result of our theory is somewhat 
worse than that of the RPA, although we have included 
additional spin correlations. This behavior is analogous 
to the results for the magnetization in Ref. [13, where 
the RPA is found to be closer to the QMC data than 
the Callen decoupling. For higher fields {h > 0.3) this 
tendency is reversed, i.e., our magnetization curves lie 
slightly below the RPA curves and closer to the QMC 
data. 

Figure |2 shows the spin polarization of a i^ = 1 quan- 
tum Hall ferromagnet measured by magnetoabsorption 
spectrosocopy^ in comparison with our theory and QMC 
data, where h — 0.32 is takeni2i4 The very good agree- 
ment gives a justification for the use of an effective 2D 
Heisenberg model to describe this itinerant ferromagnet. 
This should be further confirmed by the comparison of 
other thermodynamic quantities (magnetic susceptibility, 




FIG. 1: Magnetization of the ID (a) and 2D (b) Heisenberg 
ferromagnet in magnetic fields of strengtlis li=f .0, 0.8, 0.6, 
0.4, and 0.1, from top to bottom, as obtained by the Green's- 
function tlieory (solid), the ED (o), and the Bethe-ansatz 
method (■), compared with RPA results (dotted) and QMC 
data (•, Ref. ff). The inset shows the low-field magnetization 
of the ID model at h = 0.05 and 0.005 from top to bottom. 



specific heat) with experimental data which, however, 
not yet available. 
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B. Magnetic susceptibility 

Let us consider the isothermal susceptibility Xs — 
d{S'')/dh shown in Figs. 01 and |3| For /i = 0, Xs di- 
verges at T = indicating the ferromagnetic phase tran- 
sition. In one dimension (see inset of Fig.|3^), the Bethe- 




FIG. 2: Magnetization for the v = 1 quantum Hall ferromag- 
net calculated at h=0.32 (solid) in comparison with QMC (•, 
Ref. U and experimental data (D, Ref. [J). 



ansatz result lim Y„r 



^(Ref-ll 



is reproduced by 

the spin-rotation-invariant Green's-function theoryiSS At 
low temperatures, the deviation of the ED data, calcu- 
lated from Xs ~ r~^ / ^ C^, is ascribed to finite-size 

R 

effects. Contrary, the RPA curve, obtained by Eq. I|25|) 
which yields A = Vl + 4T2 - 1 (Ref.ll^), strongly devi- 
ates from the exact result. In the 2D system, the zero- 
field susceptibility (see inset of Fig. 01 of the rotation- 
invariant theory corrects the numerical values given in 
Ref. 20 and agrees well with the ED result. At non-zero 
fields we have XsC^ = 0) = 0. Therefore, the suscep- 
tibility has a maximum at T^, where T^ increases and 
Xs{T^) decreases with increasing field. In one dimension 
(Fig. O, the good agreement between Green's-function 
theory, Bethe-ansatz method, and ED corresponds to the 
results depicted in Fig. QJi. The deviation of our theory 
for the 2D model at /i = 0.4 and T < 1 from the ED data 
(Fig.^ is due to finite-size effects in (S*^), as can be seen 
in Fig.^. 

Recently, the field dependence of the position of the 
susceptibility maximum has been discussed in connec- 
tion with experiments on Lao.91Mno.95O3 showing a shift 
of the maximum in the temperature derivative of the 
electrical resistivity at an applied field according to h^^^ 
(Ref. \2m . Assuming that this maximum coincides with 
the maximum in the susceptibility due to spin scatter- 
ing, the dependence T^{h) was investigated in terms of 
Landau's theory, developed for anisotropic systems with 
Tc{h = 0) ^ 0, which yields T^ ex h^'^ (Refs-IHandlll. 
Within Landau's theory Sznajd— claims that this power 
law also holds for isotropic ferromagnets in a field. Con- 
sidering, as a further characteristics, the height of the sus- 
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FIG. 3: Isothermal susceptibility of the ID ferromagnet at 
low fields (a), h = 0.005 and 0.05, from top to bottom, and 
at higher fields (b), h=0.4, 1.0, and 2.0, from top to bottom, 
where the Green's-function (solid), ED (o), Bethe-ansatz (■), 
and the RPA results (dotted) are shown. In the inset the ID 
zero-field susceptibility is depicted. 



ceptibility maximum XsC^m)' Landau's theorjii^ yields 
Xs(T^) ex m-^{T^) oc /i'^/^. In Ref. [H the isotropic 
spin chain was investigated by a real-space renormaliza- 
tion group method and T^ ex h'^ with 7 = 0.696 for 
0.1 < h < 5 was found; however, XsC^m) ^^^ ^^ot calcu- 
lated. 

To analyze the power-law behavior in more detail, the 
dependence T^i{h) is plotted logarithmically in Fig. [5^ 
for h ^ 0.1. Both the results of our theory and the ED 
data for T^ may be well fit to power laws in the ID (2D) 
model for h > 0.2 (0.6). The Green's-function theory 



FIG. 4: Isothermal susceptibility of the 2D ferromagnet at 
h=0.4, 1.0, and 2.0, from top to bottom. The Green's- 
function theory (solid) is compared with ED data (o) and 
RPA results (dotted). The inset shows the 2D inverse zero- 
field susceptibility. 



yields 



with 



T,l =ah^ 



1.486 



and 7 



0.739 ; ID 
0.620 ; 2D 



(26) 



(27) 



The ED results for T^ (for clarity, the fit of the ED 
data is not drawn) deviate only slightly from Eq. H26|l : 
in the ID (2D) model we obtain a = 1.051 (1.460) and 
7 = 0.765 (0.643). Considering the power-law behavior 
in the low-field region, the numerically most reliable data 
are provided by the Bethe-ansatz solution. For the ID 
system at ft, ^ 0.1, the Bethe-ansatz results are described 
by Eq. ^ with 



0.765 and 7 = 0.576 



(28) 



Comparing our results for the ID model with Ref. lla 
(7 ~ 0.7), the 7 values are in rough agreement, whereas 
the absolute values of T^ found in Ref. [H exceed our 
results by a factor of about 2.5. The dependence of the 
maximum position on the dimensionality could be used 
for the interpretation of experimental data. 

Our results for the maximum value Xs O^m) ^^ function 
of h in the same field region as before are plotted in 
Fig-inb- Again, they may be described by power laws for 
h > 0.2 (0.6) in the ID (2D) model. From the Green's- 
function theory we obtain 



xjm = bhP 



(29) 




FIG. 5: Logarithmic plot of the field dependence of the po- 
sition (a) and height (b) of the susceptibility maximum ob- 
tained by the Green's-function theory (a) and fit to a linear 
dependence (solid) in comparison with ED (o) and RPA re- 
sults (D). 



with 



0.202 
0.185 



and (3 



-0.951 
-0.914 



ID 
2D 



(30) 



The ED results for the ID (2D) model are given by 5 = 
0.206 (0.191) and /3 = -0.964 (-0.935). At low fields 
we consider, as above, only the Bethe-ansatz solution for 
the ID model which, at /i ^ 0.1, yields Eq. ^^ with 



b = 0.208 and /3 = -0.952 



(31) 



Note the remarkable agreement of the Bethe-ansatz re- 
sults for the field dependence of the maximum height 
with the findings of the theory and the ED data. 



The results obtained for the exponent /3 strongly devi- 
ate from the /3 value of Landau's theory, /3 = — | . More- 
over, we get m(T^) -const. {{S^)(T^) ~ 0.3) which 
also contradicts the law m{T^) (x h}'^ . This reflects the 
fact that Landau's theory does not hold for ID and 2D 
isotropic ferromagnets, but is valid only on the assump- 
tion of a finite critical temperature at /i = which, how- 
ever, is not realized in the ID and 2D systems. Therefore, 
the approximate agreement of the obtained 7 exponents 
with 7=1 seems to be accidental. 



C. Specific heat 

First let us consider the NN spin correlation functions 
entering the internal energy [cf. Eq. H23|l ] . which are de- 
picted for the ID and 2D cases in Figs. Eland [7| respec- 
tively. In the ID model we obtain a very good agreement 
with the ED data. On the contrary, the RPA results are 
unsatisfactory; in particular, the longitudinal correlators 
at low fields and temperatures, obtained from the exact 
representation of the internal energy, Eq. (|23|) , are neg- 
ative being incompatible with theferromagnetic SRO. In 
the 2D model at low fields, we ascribe the deviations of 
our analytical curves from the ED data to finite-size ef- 
fects; in this respect, Cfg may be considered in analogy 

to (5-) (Fig.nt). 

Figure IHl displays the specific heat C — du/dT for the 
ID ferromagnet. At ft. = 0, the temperature dependence 
of the specific heat exhibits a broad maximum, where 
the value of the maximum position resulting from the 
Green's-function theory, T^{h = 0) ~ 0.45, agrees rea- 
sonably well with the exact result T^{h = 0) ~ 0.35 
obtained by the Bethe-ansatz and ED methods. Com- 
paring our ED data at ft = with those of Ref. [2J agree- 
ing with the Bethe-ansatz results, the additional weak 
maximum at T ~ 0.1 has to be ascribed to finite-size 
effects. The broad maximum and the strong decrease of 
the zero-field specific heat at low temperatures is qual- 
itatively reproduced b y t he Green's-function theory, as 
already shown in Ref. [22. At very low magnetic fields, 
besides the high-temperature maximum, a second max- 
imum at low temperatures develops which has not been 
reported before. The occurence of two maxima in the 
specific heat is indicated by our theory, however, at too 
high fields (0.03 < ft < 0.07). In a detailed Bethe-ansatz 
analysis, two maxima in the specific heat are found in 
the field region < ft < 0.008 (see inset of Fig. ^). 
At ft ^ 0.008, only one maximum appears. Considering 
the low-temperature maximum at ft ^ 0.01, the position 
and height C{T^ i) are given by the power laws 



rpC 



T^i = 0.596 h° 



C(Tg^) = 0.513 ft° 



(32) 



Note that the exponent of T^ ^ nearly agrees with that 
of Tx given by Eqs. ^ and' (EHl- From the low-field 
specific heat it becomes evident again that our theory 
provides an improved description of SRO, as compared 
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FIG. 6; Nearest-neighbor transverse (a) and longitudinal (b) 
spin correlation functions of the ID ferromagnet at the fields 
h=0.1, 0.4, 1.0, and 2.0, from left to right. The Green's- 
function theory (solid) is compared with ED (o) and RPA 
results (dotted). 



FIG. 7: Nearest-neighbor transverse (a) and longitudinal (b) 
spin correlation functions of the 2D ferromagnet at the fields 
h=0.1, 0.4, 1.0, and 2.0, from left to right. The Green's- 
function theory (solid) is compared with ED (o) and RPA 
results (dotted). 



with RPA (cf. Fig. |S|d) which does not yield a double 
maximuiii. 

The specific heat for the 2D model is plotted in Fig.|51 
We get a good agreement with the ED results, in par- 
ticular, as the position and height of the maximum is 
concerned. Note that the small low-temperature bump 
in the ED data for /i = 0.1 is a finite-size effect. The 
RPA curves at low fields exhibit a too large maximum 
height which we ascribe to a poor description of SRO in 
RPA (see also Fig.[71l. 

In the ID and 2D systems at the fields h > 0.4 and 
h > 0.1, respectively, the position of the specific-heat 
maximum obtained by the Green's-function theory may 



be described by the linear law 
with 



0.433 J L_ J 0.310 ; ID 
" ^ ^ 0.463 '^'^^ ° - 1 0.685 ; 2D 



(33) 



(34) 



IV. SUMMARY 

In this paper we developed a Green's-function theory 
of the ID and 2D S = 1/2 Heisenberg ferromagnet in a 
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and an 4 X 4 square lattice and obtained exact Bethe- 
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FIG. 8: Specific heat of tfie ID ferromagnet obtained by tlie 
Green's- function (solid) and Bethe-ansatz method (■) at low 
fields (a), /i = 0, 0.005, and 0.03, from bottom to top, with 
the ED data denoted by ©, o, and (g), respectively, and at 
higher fields (b), h = 0.1, 1.0, and 2.0, from left to right, 
in comparison with ED (o) and RPA results (dotted). For 
clarity, the Bethe-ansatz data for h = 0.005 and 0.03 at low 
temperatures are joined by dotted lines. For low fields the 
RPA data are not drawn because of the too high maximum 
(cf. (b)). The inset exhibits the Bethe-ansatz results for very 
low fields, /i = to 0.01 in steps of 0.001, from bottom to top. 
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FIG. 9; Specific heat of the 2D ferromagnet at at ft = 0.1, 
1.0, and 2.0, from left to right, showing the Green's-function 
(solid), ED (o), and the RPA resuhs (dotted). At ft = the 
Green's-function theory (dashed) is compared with ED data 

(0). 



ansatz results for the Heisenberg chain from an eigen- 
value analysis of the quantum transfer matrix. Stimu- 
lated by recent disscussions vi^e analyzed the field depen- 
dence of the maximum in the temperature dependence of 
the isothermal magnetic susceptibility. We found power 
laws for the position and height of the susceptibility max- 
imum which are shown not to be related to the predic- 
tions of Landau's theory. Paying particular attention to 
the specific heat of the Heisenberg chain, in a detailed 
Bethe-ansatz analysis the existence of two maxima in 
the temperature dependence of the specific heat at very 
low magnetic fields was proven for the first time. The 
field dependences of the position and height of the low- 
temperature maximum obey power laws. Altogether, we 
analyzed the effects of dimensionality (ID vs 2D) on all 
thermodynamic quantities which may be relevant for the 
comparison with experiments. 
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